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Abstract: Using a four dimensional approach, we show that the singularities for 
small gluon momenta, which arise in the usual three dimensional treatment of the anni- 
hilation decay, disappear if all poles in the relative energy are taken into account correctly 
in the integration. We obtain an explicit formula for the decay width which involves a 
non locality originating from the kinetic energy. We calculate not only the familiar log- 
arithmic dependence on the binding energy, but also the constant to be added to the 
logarithm. For positronium this differs from an earlier result in the literature and leads 
to a modified life time. In QCD there is a non abelian loop graph which contributes to 
the decay amplitude to the same order as the tree graph. Contrary to the original ver- 
sion, which contained some errors, the constant turns out to be large and negative. For 
reasonable parameters, this negative constant practically cancels the positive logarithm, 
so that the whole approach of expanding with respect to the binding energy breaks down. 
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1 Introduction 

It is known for a long time that the usual non relativistic approach for calculating the 
decay width of the 1 P 1 +_ -state into three gluons, as well as the annihilation of 3 P 1 f + into 
a gluon and light quark- ant iquark pairs, has infrared divergences originating from soft 
gluon momenta. In the original paper by Barbieri, Gatto, and Remiddi pQ the problem 
was overcome by introducing a fictituous width for the decay into a real and a virtual 
photon (or gluon). Their result was (i?2i( r ) is the radial wave function) 



Actually only the singular logarithmic dependence on the binding energy could thus 
be calculated, while the constant above was simply ignored. This is a shortcoming 
because such a constant, even if not particularly large, could be quite important since 
the logarithm is not large in QCD and even only about 12.6 in QED. 

In a series of papers Bodwin, Braaten, Lepage [2] and collaborators clarified some 
aspects of the problem by applying the methods of non relativistic QCD (NRQCD) 3J. 
The quarkonium state is not only made up of quark and antiquark (QQ) but contains also 
Fock state components with quark-antiquark-glue (QQg) etc. Although this admixture of 
glue is suppressed with respect to the dominant QQ-sta.te, it contributes to the decay with 
the same order. The reason is that the quark-antiquark annihilation in the QQ-st&te is 
also suppressed, because the QQ wave function at the origin vanishes in the P-state. The 
QQ in the QQ g-componet, on the other side, can be in an S-state, thus compensating the 
suppression of the Fock state component. The divergences found previously are canceled 
and the part with the logaritmic dependence on the binding energy can be calculated, 
the constant which accompanies the logarithm remains, however, unknown. 

In a further paper on the subject the present author jl] calculated the glue content in 
heavy quarkonia in order to obtain more definite statements. Although this calculation 
was succesful and gave reasonable results, it turned out that it was, at least up to now, 
not helpful for calculating the decay width. 

In the present work we go back to a more fundamental formalism. It involves the 
four dimensional T-matrix for the annihilation of the quark-antiquark pair into gluons, 
as well as the four dimensional Bethe Salpeter wave function, taken in the approximation 
coming from the static Salpeter equation. Both quantities depend on the relative energy 
P° = (Pi ~P2)°/2 of the quark-antiquark pair. The usual non relativistic approach is 
equivalent to considering only the poles of the BS wave function when integrating over 
p°. The p°-dependence of the T-matrix is, however, ignored by fixing the quarks on mass 
shell there. As will be discussed in the following, this procedure is legitimate in most 
cases, like decays of S-wave states, or decays into two real gluons. It fails, however, for 
the decay of the 1 P 1 H -state into three gluons if one of the gluons is soft. 

The contributions from the poles in the T-matrix are only of relevance if one of 
the gluons is soft. In this region, however, they are crucial, because they cancel the 
singularities in the leading term. Thus the correct consideration of all poles gives an 
expression which is free of divergences. Another phenomenon which was already found 
in |3] arises in the soft region: A non abelian graph, where a soft transversal gluon splits 
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into two Coulomb gluons, while the latter couple to quark and antiquark, contributes 
to the same order as the tree graph. The reason is that the direct coupling of a soft 
transversal gluon to a non relativistic quark is also suppressed. 

Our result is not directly expressible by the derivative of the radial wave function at 
the origin, because it contains expressions of the form k + \ E\ +p 2 /m in the denominator, 
with k a gluon energy, E the binding energy, and p 2 /m the kinetic energy. The presence 
of the kinetic energy introduces a non locality. For a Coulomb potential one can carry 
through the calculation analytically till the end. We agree with the formulae in the 
literature as far as the logarithmic dependence on the binding energy is concerned. 
Beyond this we can calculate the constant term to be added to the logarithm. For the 
decay of the singlet P-state of positronium into three photons (where the non abelian 
contribution is of course absent) our constant differs from the one given by Alekseev 
[3]. In QCD the resulting constant is large and negative for reasonable parameters, and 
practically outweighs the positive logarithm. This means that the whole approach breaks 
down unless ot s is very small. 

In our approach no QQg component in the wave function shows up at all. All 
contributions come from the four dimensional QQ wave function alone. Nevertheless 
one may reinterprete our result in terms of the non relativistic picture and the presence 
of a quark- antiquark glue component in the wave function. 

In sect. 2 we discuss the situation and derive the general formalism. In sect. 3 we 
specialize to a Coulomb wave function and present an analytic formula. In sect. 4 we 
discuss the experimental status and future applications. 



2 Four dimensional versus three dimensional approach 

Consider an annihilation decay of a quarkonium state in the rest frame into three (or 
analogous into two) gluons. In the familiar three dimensional approach one proceeds as 
follows. The quark energies p\ and p\ i n the T-matrix for QQ — > gluons are taken on 
mass shell and the quark momenta smeared with the Schrodinger wave function. The 
amplitude, with the gluon momenta defined as in Fig. 1, then becomes 

< mK >= + jo f T(pm ^ (21) 

^/(27r) 3 m J 

with p^ = (pi — p2) fJ '/2 the relative momentum. The dependence of the T-matrix on 
the gluon momenta ki, k2, k3, polarization vectors e±, 62, 63, and color indices a, b, c has 
been suppressed. We used the standard covariant normalization of states and relativistic 
normalization of T, the Schrodinger wave function is normalized to 1. The T-matrix for 
quark- antiquark into gluons is a slowly varying function of p, with a scale set by the 
quark mass m, while the wave function ip(p) is non relativistic, i.e. dominated by soft 
momenta of the order of a s m. Therefore one can expand the T-matrix around p = 0. 
For S-states it is sufficient to consider the leading term T(0), the remaining integral gives 
the wave function in position space at the origin, 

T(0) J 4>(p)d 3 p = T(0)(2n f 2 m- (2-2) 
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For P-states (with magnetic quantum number m) one has to expand up to first order in 
p, resulting in 



[(Vp)"T( P )] 

-*[(v P ) 



p=0 



p n ip(p)d 3 p = 

Y lm (e (n) )(27if 2 R'(0), 



p=0 



(2.3) 



with R(r) the radial wave function. 

For decays into two gluons, which give the leading contribution for states with pos- 
itive charge conjugation, i.e. 1 S (j ~ + or 3 Pj + for J = 0, 2 ( 3 P^ + cannot decay into two 
real gluons) there are no problems. Both gluons have to be hard from kinematical rea- 
sons. Therefore the denominator of the quark propagator in T(p) cannot vanish and no 
singularities can arise. 

Let us next discuss the three gluon decay shown in Fig. 1 which is the leading 
contribution for annihilation decays of states with negative charge conjugation, i.e. 3 S{ ~ 
or 1 P 1 +-. 




Figure 1: The lowest order contribution to the three gluon decay. 

With the quark energies taken on mass shell, and the quark momenta considered only 
up to first order, the denominators of the two quark propagators become (we always use 
kj = |kj| in the following) 

(pi — ki) 2 — m 2 + ie = — 2mki + 2p 1 k 1 + ie, 

(k 2 - P2) 2 -m 2 + ie = -2mk 2 + 2p 2 k 2 + ie. (2.4) 

This is the well known reason for the possible appearence of divergences which can arise 
from small ki or k 2 . 

Let us now view the situation in the framework of a four dimensional treatment of 
the three gluon decay, focusing on the problematics for the P-state. The amplitude in 
Fig. 1 has the form 

< KkiKlK >= wty *.+*»-*)„ (2 5) 

'(2-Kfm 



J Tr[T(p,p°)T(p) 



(\E\ + p 2 /m) ijj(p) dp d 3 p 



(p° + \E\/2 + p 2 /2m - ie) (jpP - \E\/2 - p 2 /2m + ie) 
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As before, T denotes the T-matrix for the annihilation of the quark- ant iquark pair into 
gluons. T(p) is the spin wave function of the bound state. For a singlet state it reads 
T(p) = A(pi)(7 5 /v / 2)A c (p 2 ). Here A(p fc ) are the projectors to positive energy states, in 
our case we need them only up to order p k where they read 

Finally is again the Schrodinger wave function, while E is the binding energy. We 
used the well known p°-dependence of the Bethe-Salpeter wave function as it arises from 
the static approximation to the BS equation. 

We first discuss the trace in the numerator. For the graph in Fig. 1 it is 

Tr{r(p) 7 " 2 [^ 2 - & + m] 7 ns [(&- A) + m] 7 " 1 }. (2.7) 

This is only needed in zeroth order of p for S-states, and up to first order for P-states. The 
zeroth order of the trace is proportional to kik 2 (if ultrasoft terms are dropped), therefore 
there is no trouble for the three gluon decay of S-states; the factors of the numerator 
cancel the singularities in ki and k 2 of the denominator. This is no longer the case for 
the terms of order p. The trace up to order p, contracted with the polarization vectors, 
is 

2i{ ((A* + \E\/2 - p°) (k 2 + \E\/2 + p°) + (ki • k 2 )) [ei x e 2 ] • e 3 

+ (ei • e 2 )[ki x k 2 ] • e 3 + (ei • k 2 )[e 3 x e 2 ] • k x - (e 2 • ki)[e 3 x ei] • k 2 } 

+4i{(e 2 • p)[ei x e 3 ] • k x + (ei • p)[e 2 x e 3 ] • k 2 }. (2.8) 

We next have to perform the integration over p° in (2.5). The usual approach is 
equivalent to fixing the energies p\ and p 2 on mass shell in the T-matrix, which means 
putting p° = \E\/2 + p 2 /(2m) in the first, and p° = —\E\/2 — p 2 /(2m) in the second 
quark propagator. Only the p° dependence from the Salpeter wave function is considered 
and the p°-integration performed, thus ending up with the simple formula (2.1). In this 
way one has made three approximations: 

• The values of p° in the T-matrix have been fixed by the on shell prescription instead 
of using the residue at the pole of the wave function. 

• The contributions from the remaining poles have been neglected. 

• For P-states, where one has to expand up to order p, the p-dependence of the 
projection operators A has been ignored. 

Let us come to the correct p°-integration in the four dimensional formula (2.5). The 
Salpeter wave function has a pole e+ in the upper half plane, and a pole e~ in the lower 
half plane. There are four more poles from the two propagators in T, which we denote 
by ef , and ef. Altogether the poles are 



4 



4 = T(\E\/2 + p 2 /(2m)-ze), 

ef = —m + \E\/2 + ki=p (yjm 2 + (k x - p) 2 - ie ) , (2.9) 

ef = m - \E\/2 -k 2 T (\/™ 2 + (k 2 + p) 2 - «e ). 

We have to calculate the following integral where the first four factors come from 
the quark propagators in the T-matrix, the last two from the four dimensional wave 
function: 



_J_ r°° (\E\ +p 2 /m) dp 

2m J-oo (p« - e+)(p° _ e -)( p o _ e + )(p o _ er)(p o _ e + )(p o _ e - } ■ K - ) 

We perform the p°-integration by closing the integration contour in the upper half plane 
(which, temporarily, introduces an apparent asymmetry between quark and antiquark). 
There are three contributions originating from the three poles at e+, ef, e\: 

, = -(| B | + P 7 m )(J- + i- + i-). (2.11) 

The denominators which appear here are 

D w = li(e+-e n ), D 1 = ]X(ei-e n ), D 2 = U'(e+ - e„), (2.12) 

n n n 

where the products fin' run over the five remaining poles. 

If all three gluons are hard, the contribution 1/D W from the pole e+ is dominant 
because D w contains the small ultrasoft difference — e w = —(\E\ + p 2 /m). This 
cancels the corresponding factor in front of (2.11). All the other differences appearing 
in the three denominators are hard for hard gluons, therefore the contributions 1/D 1 
and I/D2 are suppressed compared to the leading term 1/D W by a relative factor of the 
order (\E\/m + p 2 /m 2 ). The momentum p 2 is on mass shell on the pole e^J, while p\ is 
slightly off shell. One may, however, also put p\ on mass shell in D w , which corresponds 
to replacing e+ — by — ef. This again only causes a relative error of order \E\jm 
and leads back to the old formula from the three dimensional treatment. 

The situation is more subtle if one of the gluons is soft or ultrasoft. The other two 
are then, of course, hard. 

We start with the case where k 3 — > 0, i.e. when the quark propagators in the T- 
matrix become identical and create a double pole. Both D\ and D 2 become singular for 
k 3 — > which is due to the factor 



-h - ^m 2 + (kx - p) 2 + ^/m 2 + (k 2 + p) 2 (2.13) 

(we used k\ + k 2 + k 3 = 2m — \E\). This factor enters, however, with opposite signs in 
D\ and D 2 , while all the other factors become identical for k 3 — > 0. Therefore there is 
no singularity for k 3 — > in the sum. 
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For ki — > the first propagator gets near to the mass shell. This behavior shows up 
only in D w and is due to the difference 



e+-e l = -\E\ - p 2 /(2m) +m - k x - ^Jm 2 + (ki - p) 2 

-> -(fci + |£| +p 2 /^- (ki -p)/m) for ki -> 0. (2.14) 

Finally we look at the limit k 2 — ■> 0, where the second quark propagator becomes on 
shell. In D w there is a difference which becomes small, namely 



e+ - e+ = -p 2 /(2m) - m + k 2 + ^m 2 + (k 2 + p) 2 

-> A; 2 + (k 2 • p)/m for k 2 -> 0. (2.15) 

In D 2 there are two differences which become small, namely e 2 — e+ — > — (/c 2 + (k 2 -p)/m), 
i.e. the negative of the term above, and 



e 2 e v 



-\E\ - p 2 /(2m) +m-k 2 - ^m 2 + (k 2 + p) 2 

-(k 2 + \E\ + p 2 /m + (k 2 • p)/m) for k 2 -> 0. (2.16) 



In the sum the two poles at fc 2 + (k 2 • p) jm = cancel, while the ultrasoft term 
— (k 2 + \E\ + p 2 /m + (k 2 • p)/m) remains. 

Summarizing we found the following result: The contributions of 1/Di and 1/D 2 are 
negligible if all gluons are hard. They are also negligible for k3 — >• because the poles 
cancel in the sum, as well as for \ti — > 0, because they are finite in the latter limit. For 
k 2 — > finally, the poles at k 2 = in l/D w and 1/D 2 cancel, while the pole (2.16) in 
1/D 2 survives. In total one thus obtains 

1 1 1 l + Q(|ff|/m + p 2 /m 2 ) 
D w + Dl + D 2 - D ' {2A7} 

where D is obtained from D w by replacing e+ — e 2 by — e 2 . We may also replace 
e w ~ e 2 by e~ — e 2 because this difference is always hard. This restores the manifest 
symmetry between quark and antiquark. 

The spatial momentum p is soft in the bound state described by the wave function 
■0(p)- Therefore one can put p = (as well as E = 0) in all terms where p-dependent 
terms are added to hard terms. In the two critical factors, i.e. in those which become 
small for small k x or k 2 respectively, one can use the limits in the second lines of (2.14) 
and (2.16). The terms (ki • p)/m and (k 2 • p)/m are treated in first order. This finally 
leads to the following simple form for the integral / in (2.10): 

1 1 + m m - 7% ifeix */ s + 0(\E\/m + p 2 /m 2 ) 

j _ m(ki + \E\+p z /m) m(k2 + \E\+p^ /m) VI 1/ r I J , , 

''"4m 2 " ""^^ (ki + \E\ +p 2 /m) {k 2 + \E\ +p 2 /m) ' l ' ' 

The apparent singularities for ki — > and k 2 — > 0, which showed up in the three 
dimensional approach have turned out to be spurious. They disappear if all the poles 
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in the relative energy p° are considered in the p°-integration. For j = 1,2 there remain, 
however, denominators of the form kj + \E\ + p 2 /m which become small for small kj 
and enforce to keep a finite binding energy in order to avoid a divergence. Because the 
binding energy E and the kinetic energy p 2 /m are of the same order, one must, however, 
also keep the term p 2 /m for consistency. This leads to a non local dependence upon the 
wave function. 

We have to multiply (2.18) with the trace (2.8) and the remaining contributions 
from the expression (2.5). If terms contain e.g. a k™ in the numerator, one can replace 
k% + \E\ + p 2 /m — > k± in the denominators of (2.18). This finally leads to 

/ Jabc v 

< KKK\ K >tree= ( J l^fb^^ + k 2 + k 3 - K) X 

9 3 



4VN 



7r 3/2 m 3/2 

with 



J F(p) , 0(p) d 3 p + permutations, (2.19) 



( , [ei x e 3 ] -ki(e 2 -p) 1 / (ki ■ p) _ (k 2 ■ p) \ 

lPJ fc!(fc 2 + |£| +p 2 /m) + Amk x k 2 \ k x k 2 ) 

(k x k 2 + (kt • k 2 ))[ei x e 2 ] • e 3 + (e x • e 2 )[ki x k 2 ] ■ e 3 + 2(e 1 ■ k 2 )[e 3 x e 2 ] • k x . 

The permutations are the remaining five permutations of the three gluons. The factor 
(d abc / (4\/~N)) in front is the color factor, with N = 3 the number of colors. For the decay 
of the singlet P-state of positronium into three gammas it has to be replaced by 1, while 
g is replaced by e. The first term in E(p) originates from the product of the order p term 
of the trace with the order 1 term of the denominator, the remaining contributions come 
from the order 1 term of the trace, multiplied with the order p term of the denominator. 

In QED (2.19) would in fact be the leading contribution to the decay amplitude. In 
QCD we encounter, however, the same phenomenon which we found already in |-fj when 
calculating the admixture of the quark- ant iquark-glue component to the wave function. 
For the case that k 2 , say, is ultrasoft, the loop graph in Fig. 2, where the dashed 
lines denote Coulomb gluons, contributes to the same order as the tree graph in Fig. 1 
(although it will not lead to a logarithmic dependence upon the binding energy). 

The reason is that the tree graph is suppressed by the coupling of a transverse gluon 
to a non relativistic quark. The loop graph, on the other hand, is suppressed by two 
more coupling constants, but this is compensated because the coupling to the quarks 
via Coulomb gluons is not suppressed. Other loop graphs are not relevant, because they 
don't have a sufficient number of denominators which can become small simultaneously. 
The same is true for higher order corrections in the wave function, e.g. exchange of a 
transverse gluon. Also the graph in Fig. 2 is only relevant in the region of ultrasoft k 2 . 
This simplifies the calculation considerably. 

If one first integrates over q° in Fig. 2 one finds a similar situation as before. Spurious 
poles in k\ and fc 3 which show up at some residues cancel in the sum. Only the variable 
k 2 appears in a denominator which becomes small but finite for ultrasoft k 2 . After having 
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Figure 2: T/ie /oop graph for the three gluon decay. Note that the indices of the gluons 
have been chosen such that the dominant contribution comes from k 2 soft. 



performed the g°-integration, the whole p°-dependence is in the BS wave function and 
the p°-integration can be immediately performed. Finally one can do the integration over 
d 3 q which is simplified by the fact that one can drop the k 2 in the Coulomb propagators. 
The resulting expression is 



with 



/ Jabc 

< Kk^K > loop = -(^=) li^Y^ih + k 2 + k 3 -K)x (2.20) 



g 5 N[ei x e 3 ] • k x f (e 2 ■ p) f rfm(k 2 + \E\, ^ 



\§TT^I 2 ^[m J p 3 V p 
+ permutations, 



7T X 

fix) = — arctanx. (2.21) 

M J 2 1+x 2 K ' 



3 Three gluon (or 3 7) annihilation of the 1 Pf state 

Due to the non local dependence on the momentum p in (2.19), (2.20) one cannot simply 
express the result by the derivative of the wave function at the origin. Instead one has 
to perform the p-integration explicitly. We carry this through for the case of the lowest 
Coulomb P-state wave function which reads (we drop an irrelevant factor i) 



* (p) = — ^^ ( °) (| P pl1|fl|)' y '" (W - (3 - 1) 

with \E\ = \E 2 \ = (C F a s ) 2 m/16 and C F = (N 2 - 1)/(2N) = 4/3. To make contact 
with the literature we have split off the derivative of the wave function at the origin, 
#21 (0) = 2(m|£|) 5 / 4 / v / 3. The p-integration in (2.19) and (2.20) can now be performed 
and we obtain 
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< kXKlK >c, 



/ Jabc \ 

olllomh = J t(2Tr) 4 5^(k 1 + k 2 + h- K)R' 21 (0) x 

[g 3 ftree(ki,k 2 ) + g 5 fioo P (h,k 2 )}+ permutations. (3.2) 



The functions /t ree (&i, k 2 ) and fi oop (ki, k 2 ) arise from the tree graph and the loop graph, 
respectively and read 



ftree{k\i k 2 ) — 

2V2 



+ 



[ei x e 3 ] • k!y lm (e 2 ) (8E 2 + 12\E\k 2 + 3k 2 2 - &y/\E\(\E\ + k 2 f /2 ) 



Yi m ( — ) ~~ ^1111(7- ) 



^/2m^l 2 k x k 2 



k± 



(k x k 2 + (ki • k 2 ))[ei x e 2 ]e 3 + (e x • e 2 )[ki x k 2 ]e 3 + 2(ei • k 2 )[e 3 x e 2 ] • k x 



fioop(ki, k 2 ) — 



2y/2- 



[ei x e 3 ] • kiFi m (e 2 ) 



3 • 



(3.3) 



In the calculation of the width we use the gluon energies k\ and k 2 as independent 
variables in the phase space integral, the region of integration is given by the triangle 



< h,k 2 < m- \E\j2 



m, 



k\ + k 2 > m — \E\/2 m m. 



(3.4) 



As indicated, one can neglect the ultrasoft binding energy in the boundaries. 

One next has to square the T-matrix in (3.2), perform the summation over gluon 
polarizations, and integrate over ki and k 2 . Of course one can omit the permutations in 
one factor and, instead, multiply by a factor 6. The integration over k\ gives a factor 
k 2 from the boundaries of phase space. Therefore only those terms are sensitive to the 
binding energy E where both factors become small for small k 2 (see the structure of 
(2.19)). For the tree term this is only the case in two permutations, namely the identity 
and the exchange k\ <-> k 3 , in the remainig four one may drop E and perform both 
integrations analytically. In the two potentially "dangerous" terms, one has to keep 
the binding energy. The /^-integration can be performed analytically (with the help of 
Mathematica). The remaining /^-dependent function can then be split in two terms. 
The first term can be integrated analytically and contains the logarithmic dependence 
on E, while in the second term one can drop E and subsequently also perform the k\- 
integration analytically. In the square of the non abelian term, as well as in the mixed 
term, only the identity and the permutation k\ <-> k 3 contribute. One may put E = 
there, no logarithmic singularity arises from these terms. 

For positronium, where the non abelian contribution is absent, the result for the 
decay width into three gammas finally becomes 
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r(ip+- _> 37) 



8a3 i/?' 



m ,25 n 2 „ , „ s ,3?H 3 N ,49tH 5, 



In —-(— + — + 21n2) - (— - -) + 



am 
96vr 



,B| 9 64 y v 32 4' v 192 2' 

, 16 ,25 7T 2 , , .3tt 2 3, ,497r 2 5 N 

In— - — + — + 21n2 + 

a 2 v 9 64 ; v 32 4 ; v 192 2' 



(3.5) 



The width for the quarkonium decay into three gluons becomes (5/18 is the color factor) 



T( 1 P 1 + - - 3 glnons) = T^^iW* 

, m ,25 vr 2 , , .37T 2 3 N ,497r 2 5, 9 9 

In— - — + — + 21n2 + - + — 

LB v 9 64 ' y 32 A J v 192 T 8 64 



(3.6) 



Let us first discuss the positronium case (3.5). Comparing with the formula of Alekseev 
[5] which reads 

r( 1 Pi +_ -> 37) A « e fc see ,; = 7^7— In (3-7) 

967r cr 

we see that the logaritmic term lnl6/a 2 coincides, but the constant differs. Instead of 
In 2 = 0.693 in Alekseev's formula we have -(f + g +2 In 2) - - |) + (^g- - §) = 
—4.47475. The life time according to our formula becomes r = 0.00594 s, which is a 
factor of 1.6 larger than the life time which is derived from (3.7). 

There is also a connection to another interesting problem of QED. Recently Manohar 
and Ruiz-Femema |6| reanalyzed the low energy photon spectrum of the three photon 
decay of orthopositronium in the framework of NRQCD. This work was motivated by 
the observation of Pestieau and Smith that the spectrum, as derived by Ore and 
Powell jS] long ago, violates Low's theorem jU]. The authors found agreement with the 
Ore Powell formula for the case that all photons are hard, but obtained deviations if one 
photon becomes soft. The low energy photon spectrum (the energies of the remaining 
two photons are integrated) does not behave like the Ore Powell spectrum but is much 
softer. 

One can try to understand the origin of such a suppression in our approach. An 
analogous calculation for the orthopositronium 3 S 1 would show that the k has been 
replaced by k + \E\ + p 2 /m in some of the denominators, i.e. for k — ► there is a 
suppression factor. This is, however, not sufficient to give the correct result. The point 
is, that for ultrasoft kj, electric or magnetic radiative transitions to the states 3 P + 2 + or 
1 Sq + , which subsequently annihilate into two photons, become relevant. To take into 
account these bound states correctly, one would need the exchange of a whole ladder. For 
this special problem the methods of non relativistic QED are certainly more appropriate 
(I thank P. Ruiz-Femema for an enlightening correspondence). 

We next discuss our result (3.6) for the quarkonium decay into three gluons. As far as 
the logarithmic term is concerned it coincides with the formula given by Barbieri, Gatto, 
and Remiddi pQ . Beyond this we were able to calculate the associated constant term to be 
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added. The five contributions arise from the square of the expansion (in first order with 
respect to p) of the numerator of the tree graph, the mixed term between numerator and 
denominator, the square of the expansion of the denominator, the mixed term between 
tree graph and loop graph (only the expansion of the numerator contributes), and the 
square of the loop graph. Numerically they are 

-4.31828 - 0.17528 + 0.01881 - 1.12500 + 0.14062 = -5.45913. 

To get a feeling for the importance of the constant term let's use \E\/m = (4a s /3) 2 /16 
and choose a s = 0.2 which gives m/\E\ = 225 and hx(m/\E\) = 5.4161. This shows that 
the negative constant outweighs the positive logarithm and formally leads to a negative 
widths! The whole approach breaks down and makes only sense if m/\E\ is much larger 
than the value above. (In the original version there were some errors which led to a very 
small constant. The conclusions thus have changed drastically!) 



Some remarks on the error are appropriate. It is a term of the order ^J\E\/m in 
the brackets of (3.5) and (3.6). A consistent calculation of this correction would be 
extremely cumbersome, considering the many places where we have dropped the binding 
energy. It is, however, possible to calculate the correction coming from the last step, 
when performing the integration in the two permutations which lead to the logarithmic 



term. It results in the addition of a term 13.1y\E\/m ~ 0.87 in the bracket of (3.6). 
The correction is thus not dramatic, in spite of the quare root behavior and the large 
factor in front. 

One might be tempted to simplify the whole procedure by replacing the non local 
p 2 /m-dependence in (2.18), (2.19) by some average p 2 /m. This would directly lead to a 
result involving |i?'(0)| 2 for any wave function. In the non abelian loop graph there is a 
dependence on the variable (p — q) 2 , and it is hard to find a reasonable approximation 
for this. Therefore we restrict the discussion to the abelian case. The result would be 



r( ip+- _» 3 7) 



a 8 mf, 16 /5 7tt 



96vr 



ln—-(--— + \n(l + P y(m\E\) 



(3.8) 



One could next use the virial theorem for the Coulomb potential and put p 2 /m = \E\. 
This would obviously lead to a constant different from that in (3.5). The width would be 
multiplied by a factor of 1.33 as compared to (3.5). Thus one obtains a bad approximation 
even in QED where the logarithmic term ln(l/a 2 ) dominates. This result shows that 
the nonlocality is essential. On the other hand it demonstrates the stability of the 
logarithmic term ~ ln(l/a 2 ) against any reasonable approximation which respects the 
fact that p 2 /m and E are of the same order. 



4 Conclusions 

To our knowledge there are no data on the decay of the singlet P positronium state. For 
the cc-system there is only a doubtful candidate [TDj at 3526 MeV for the 1 P 1 + ~ state h c , 
with width < 1.1 MeV and no annihilation decays observed. For the bb system there is 
not even a candidate. This is unfortunate because these states would not only provide 
an excellent further test for the nature of the long range spin dependent potential, but 
also for our ideas about annihilation decays. To do this one has, however, to go beyond 
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the approach of expanding with respect to the binding energy. Our prediction for the 
annihilation width is a modification of that given in pQ and reads 



rej?- ) = iOg. (ln ^_ 5 459) (4 1) 

This shows again that the whole approach can only be trusted for very small binding 
energies. 

From a purely theoretical point of view the annihilation decays under consideration 
are highly interesting and provided some puzzles in the past. We have shown here that 
a four dimensional treatment of the decay immediately leads to a result which is free of 
infrared diverences and which, furthermore, allows a definite calculation for any given 
wave function. Notwithstanding the fact that NRQCD has certainly greatly improved 
our understanding of low energy QCD and has led to a systematic and effective way 
for calculations, we believe that for some problems, like the one discussed here, a four 
dimensional treatment is more transparent and effective. 

We derived our result by using the QQ- wave function only. Nevertheless it is easy to 
make contact with the three dimensional approach. The contributions from the poles in 
the T-matrix can, alternatively, be interpreted as contributions from QQg components 
(with a soft gluon) in the wave function. The advantage in our approach is that we 
obtain this contribution immediately in an explicit form. Obviously these two ways of 
viewing the situation would also apply to decays into an arbitrary number of gluons. 

The decay of the spin triplet state 3 P 1 ++ into a real gluon plus a virtual gluon which 
further decays into a light quark-antiquark pair shows the same infrared problems as the 
decay discussed here. Phenomenologically it is more interesting, because the charmonium 
state Xci(l-P) is experimentally well established and theoretical results can be compared 
with the data. This will be undertaken in a forthcoming work. 

Acknowledgement: I thank D. Melikhov for reading the manuscript and for valu- 
able suggestions. 
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